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We present a determination of the nucleon parton distribution functions (PDFs) and of 
the strong coupling constant a 3 at next-to-next-to-leading order (NNLO) in QCD based 
on the world data for deep-inelastic scattering and the fixed-target data for the Drell-Yan 
process. The analysis is performed in the fixed-flavor number scheme for nf =3,4,5 and 
uses the MS scheme for a s and the heavy quark masses. The fit results are compared with 
other PDFs and used to compute the benchmark cross sections at hadron colliders to the 
NNLO accuracy. 

The nucleon PDFs play crucial role in the collider phenomenology and very often they put 
a limit on theoretical prediction accuracy, particularly for the calculations in the next-to 
-next-to-leading order (NNLO) in QCD. To meet 
quick accumulation of the data and steady 
progress in reduction of the systematic uncertain- 
ties in the LHC experiment we provide the NNLO 
nucleon PDF set with improved accuracy [5]. 
These PDFs are obtained from the updated ver- 
sion of the ABKM09 analysis [5] performed in the 
fixed-flavor number (FFN) scheme with the num- 
ber of fermions taken as n/ = 3,4,5, depending 
on the process used to constrain the PDFs. In the 
present analysis we replace the inclusive neutral- 
current (NC) DIS data of the HI and ZEUS experi- 
ments by the combined HERA data set, which are 
obtained from merging those of separate experi- 
ments jf . The data are substantially improved by 
cross-calibration of the separate experiments and 
by merging both statistical and systematic errors. 
Due to these improvements the combined HERA 
data provide a better constraint on the small- a; 
gluon and quark distributions. We also add to 
our analysis the inclusive charged-current (CC) 
DIS HERA data obtained by merging the HI and 
ZEUS samples. The CC HERA data provide a 
supplementary constraint on the PDFs helping to 
disentangle the small- x quark distributions. Fi- 
nally, we include the HI data obtained in a spe- 
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Figure 1: The data on F L versus x obtained 
by the HI collaboration [l] confronted with 
the 3-flavor scheme NNLO predictions based 
on the different PDFs (solid line: this analy- 
sis, dashes: JR09 12 , dots: MSTW [3]). The 
NLO predictions based on the 3-flavor NN21 
PDFs [3] are given for comparison (dashed 
dots). The value of Q 2 for the data points 
and the curves in the plot rises with x in the 
range of 1.5 -f- 45 GeV 2 . 
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cial HERA run at reduced collision energy, which are particularly sensitive to the contribution 
of longitudinal structure function Fl at small x 1. . This run was motivated by a particular 
sensitivity of the small- x Fl to the resummation effects and collinear factorization violation. 
Besides, Fl is quite sensitive to the gluon distribution therefore the data of Ref. [T] can help 
to consolidate the small- x gluon distributions provided by different groups, cf. Fig. [I] 

In our analysis the DIS data are described within the 3-flavour FFN scheme, as well as in the 
ABKM09 case. However, in the present fit we employ the heavy-quark Wilson coefficients with 
the MS definition for the c- and 6-quark masses, as suggested in Ref. [S]. For the case of MS 
definition the perturbative stability of the calculations is substantially improved. Moreover, 
in this case the constraints on the heavy-quark masses coming from the e + e~ data, which are 
commonly obtained in the MS definition, can be consistently imposed in the PDF fit. This 
leads to a reduction of the PDF uncertainties due to the heavy-quark masses. In particular, 
the errors in the 4(5)-flavour heavy-quark PDFs, which are generated from the 3-flavour ones 
using the matching conditions, are significantly improved as compared to the earlier ABKM09 
PDFs, cf. Fig. [2j 




Figure 2: The charm- (left) and the bottom-quark (right) PDFs obtained in the fit: The dotted (red) 
lines denote the ±lcr band of relative uncertainties (in percent) and the solid (red) line indicates the 
central prediction resulting from the fit with the running masses taken at the PDG values 9 . For 
comparison the shaded (grey) area represents the results of ABKM09 [6|. 

The value of strong coupling constant a s (Mz) is determined in our fit simultaneously with 
the PDFs. This approach provides a straightforward treatment of their correlation that is 
important for calculation of the uncertainties in the hadronic cross section predictions. At 
NNLO the ABM11 fit obtains the value of a s (M z ) = 0.1134±0.0011(exp.). This is comparable 
with our earlier determination a s (Mz) = 0.1135 ± 0.0014(exp.) [6], while the error is improved 
due to more accurate data employed in the present analysis. It is also in a good agreement with 
a s (Mz) = 0.11411q qq22 obtained in the analysis of the non-singlet DIS data with account of 
the QCD corrections up to the N 3 LO 10 . In the ABM11 analysis the value of a s is constrained 
both by the non-singlet and the singlet DIS data, cf. Fig. [3] For the kinematics of the SLAC 
and NMC experiments the x 2 -profile is sensitive to the power corrections including target 
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mass effects and the dynamical twist-4 terms. The latter are poorly defined by the strong 
interaction theory and therfore put limit on the accuracy of a s determined in our fit. On the 
other hand, the BCDMS and HERA data are insensitive to the power term due to kinematics 
peculiarities. Moreover, these data sets provide complementary constraints in determining 
a s [IT]. Performing the NNLO variant of our fit with the SLAC and NMC data dropped 
we obtain a s (Mz) — 0.1133 ± 0.0011(exp.), which is not affected by the power corrections. 
Furthermore, it is in nice agreement with one obtained in the nominal ABM11 fit that gives 
confidence in the consistent treatment of the power terms in our analysis (cf. also discussion in 
Ref. [12]). 
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Figure 3: The X 2 -profile versus the value of a s (Mz), for the separate data subsets, all obtained in 
variants of the present analysis with the value of a s fixed and all other parameters fitted (solid lines: 
NNLO fit, dashes: NLO fit). 

Predictions for the charged-lepton asymmetry and the inclusive jet production cross sections 
at the energy of LHC are in a good agreement with the first data collected by the CMS and 

ATLAS experiments [H CHI [TH [H] , cf ' Fi S s [UEl des P ite these data are not used in ABM11 fit. 
Moreover, the value of a s = 0.1151 ± 0.0001 (stata.) ± 0.0047(sys.) extracted from the ATLAS 
data of Ref. [TB] in the NLO [T7] is in agreement with our results. In contrast, the jet Tevatron 
data go above our predictions and the large- X gluon distribution rises significantly once they 
are included in the analysis. Note that the MSTW PDFs systematically overshoot the LHC jet 
data (cf. Fig. [5]) as well as other PDFs tuned to the Tevatron data [HI [16]. On the whole, this 
leads to the conclusion that the LHC data prefer softer gluons as compared to the Tevatron 
case. 

The Higgs production rates at the LHC and Tevatron are widely defined by the gluon 
distribution shape and the value of a s . The NNLO predictions for the cross section of Higgs 
production in the proton-proton collisions at the LHC energies calculated with different NNLO 
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Figure 4: The data on charged- lept on asymmetry versus the lepton pseudo-rapidity r\ obtained by the 
ATLAS [13] (left panel) and CMS [14] (right panel) experiments compared to the NNLO predictions 
based on the DYNNLO code [T5] and the ABM 11 NNLO PDFs with the shaded area showing the inte- 
gration uncertainties. The ABKM09 NNLO predictions are given for comparison by dashes, without 
the integration uncertainties shown. 



PDFs are displayed in Table [] At smaller collision energy, when the production rate is more 
sensitive to the large-- gluon distribution tail, the ABM11 calculations go lower than the 
MSTW08 and NN21 ones, while at high energies the difference between the predictions is 
smaller. 
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Table 1: The total NNLO cross sections in pb for Higgs production in the gluon-gluon fusion obtained 
with different PDF sets at the mass of Higgs boson Mh = 125 GeV. The errors shown are the scale 
uncertainty are based on the shifts /i = mjj/2 and fi = 2m,H and the ler PDF uncertainty, respectively. 



In summary, we have produced the new NNLO PDF set with improved accuracy at small 
x due to new input from the HERA data and refined theoretical treatment of the heavy-quark 
electro-production in the running-mass definition. The predictions based on these PDFs are in 
a good agreement with the first LHC data, which can be used in future to improve the PDF 
accuracy further. A benchmarking w.r.t. to other PDFs is performed; the differences found 
can be also reduced with the help of new HERA and LHC data. 
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Figure 5: Cross section data for 1-jet inclusive production from the ATLAS collaboration [IB] 
(left panel) and from the CMS collaboration [TS] (right panel) as a function of the jets transverse 
energy ET for hr — fip — Et compared to the result of the present analysis (solid) and 
to MSTW [3] (dashed). The theory predictions include the NNLO threshold resummation 
corrections to the jet production. 
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